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Abstract—Design and optimization (D&O) problems in 
applied electromagnetics, in particular antenna D&O, 
often rely on global search and optimization 
metaheuristics based on Nature-inspired metaphors. 
Most are stochastic in nature because the underlying 
process is, for example the swarming behavior of fish and 
birds (Particle Swarm Optimization) or the foraging 
behavior of ants (Ant Colony Optimization). While these 
algorithms generally work well (enough), their inherent 
randomness frequently limits utility in solving real-world 
antenna problems. A better approach is to use an 
inherently deterministic optimizer or to implement a 
stochastic algorithm so that it is effectively deterministic. 
This article discusses in detail Central Force 
Optimization (CFO) which is a deterministic algorithm 
analogizing gravitational kinematics. It provides an 
overview and emphasizes CFO's utility in solving the 
real-world problems routinely encountered by practicing 
engineers with examples drawn from antenna 
optimization. Because CFO is inherently deterministic 
only a single run is required, which can be a major 
advantage in addressing electromagnetic problems. Even 
in cases where the balance between decision space 
exploration and exploitation is better achieved using a 
stochastic approach, CFO offers the advantage of being 
easily hybridized to include pure or pseudo randomness. 
One method of injecting pseudo randomness is using x 
fractions which are discussed in the companion paper 
(Part II). 


Index Terms—Central Force Optimization, CFO, m 
fractions, Pseudo randomness, Global Search and 
Optimization, Metaheuristic, Antenna Design, Antenna 
Optimization, Algorithm. 


I. INTRODUCTION 


Central Force Optimization (CFO) was introduced in 2007 
[1] because the author was dissatisfied with the stochastic 
algorithms often used for antenna design or optimization 
(D&O). The term "design" refers to meeting minimum 
designer-specified performance criteria, while "optimization" 
refers to determining the "best" value of a designer-defined 
fitness function or figure-of-merit. This article focuses on 
optimization, but its concepts are equally applicable to 
design problems. 


Optimization problems fall broadly into two categories: 
benchmarks and real-world. This discussion addresses the 
real-world optimization problems that practicing engineers 
routinely encounter. A benchmark problem comprises a 
known objective (fitness) function with a known solution 
(see, for example, [67-69]), whereas a real-world problem 
requires both (i) formulating from scratch a suitable fitness 
function and then (ii) determining its unknown solution. 
Defining a good objective function can be a daunting task 
invariably compounded by a stochastic optimizer. If the 
optimizer itself generates random results there simply is no 
good way of knowing whether one objective function is 
better than another. A deterministic algorithm like CFO 
mitigates this problem by always returning the same results 
for the same setup because there is no inherent randomness. 


CFO also provides the additional advantage of allowing 
hybridization by including some measure of pure or pseudo 
randomness if doing so is beneficial. For example, purely 
deterministic CFO might be used first to decide which of 
many candidate objective functions is best followed by a 
randomized implementation to improve decision space (DS) 
exploration. Another important advantage of determinism is 
that different algorithm implementations can be compared 
quickly and with certainty. Taking CFO as an example, 
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investigating whether or not variable "gravitation" improves 
its performance is easily accomplished: only two runs are 
required, one with and another without variable gravitation. 
Stochastic algorithms, on the other hand, do not permit direct 
comparison because they require many runs to generate 
statistical performance data. And even then the questions of 
which objective function is best or which algorithm 
implementation is best cannot be answered with certainty. 


This paper is organized as follows: Section II provides an 
overview of CFO. Section III discusses generally why 
determinism is important using a typical benchmark and real- 
world examples. Sections IV and V provide real-world 
examples that illustrate the difficulties in formulating 
suitable fitness functions: a broadband HF monopole and a 
wireless power transfer (WPT) multi-arm helix. Section VI 
introduces basic CFO theory. Section VII provides CFO 
pseudocode and discusses methods for creating Initial Probe 
Distributions (IPD's). Section VII is the Conclusion. 


II. CFO OVERVIEW 
A. CFO's Gravitational Kinematics Metaphor 


CFO analogizes gravitational kinematics, the branch of 
physics that deals with the motion of masses under the 
influence of gravity. In physical space the fundamental 
equations are Newton's Law of Universal Gravitation and his 
Second Law of Motion. Imagine a group of small "probe" 
satellites flying through space, say, the Solar System. As 
these very small objects approach a massive planet, for 
example, Jupiter, they may become gravitationally trapped in 
orbit around it. Trapping generally requires a loss of energy, 
but under certain conditions, specifically "resonant returns," 
energy is conserved [52-55], and the trapped mass eventually 
escapes orbit. It seems reasonable to speculate that over time 
more and more probes would become trapped around the 
more massive planets than around the less massive ones, that 
is, the group of satellites eventually "converges" on the 
planet with the greatest gravitational pull. This is the process 
metaphorically implemented by CFO. 


CFO generalizes Newton's laws in 3D physical space to an 
optimization problem's N; -dimensional decision space [56- 


59]. The most massive object ("largest planet") is the 
optimum (maximum) value of the objective function (note 
that unlike most metaheuristics CFO performs maximization, 
not minimization). Each of CFO's "probes" is assigned a 
user-defined "mass" computed from the objective function's 
value at each sample point in the decision space. Newton's 
laws are generalized to define two equations of motion 
governing each probe's trajectory through the optimization 
problem's "landscape" (the objective function's topology over 
the decision space). As CFO progresses its probes converge 
of the objective function's maximum value, and the run 
terminates when some user-specified criteria are met. 


B. Exploitation vs. Exploration in CFO 


Because gravitational kinematics is perfectly deter- 
ministic, so too is CFO. Every CFO run with the same setup 
yields precisely the same results, this in stark contrast to 
stochastic algorithms that invariably yield different results 
one run to the next. But there is a price to be paid for 
repeatability: deterministic algorithms frequently do not 
traverse the decision space as effectively as stochastic ones 
(§2.1.1 in [60]). Determinism encourages exploitation, 
converging quickly on a solution using current data, but 
often at the expense of exploration, seeking still better 
solutions in other regions of the decision space. While 
generally it is difficult or impossible to add determinism to a 
stochastic algorithm because the underlying equations are 
fundamentally stochastic, CFO has the significant advantage 
of being readily randomized. This characteristic allows the 
algorithm designer to strike a balance between DS 
exploitation and exploration by creating hybridized versions 
of CFO. 


One effective approach to hybridization is to include 
pseudo randomness in CFO's Initial Probe Distribution [61]. 
In a truly stochastic algorithm each random variable (RV) is 
computed from a probability distribution and consequently is 
a priori unknowable. This is fundamentally different from a 
pseudo random variable (PRV) whose value is precisely 
known but arbitrarily assigned. A PRV can be specified in 
advance, for example an arbitrary sequence of numbers like 
the fractions [62] discussed in the companion article Part II 
[70], or it can be calculated in a prescribed manner. The 
PRV’s "randomness" derives from its being uncorrelated 
with the decision space topology, not that it is uncertain in 
the sense of a true RV. Because each PRV is known with 
absolute precision, either explicitly or by calculation, CFO’s 
probe trajectories remain deterministic. Including PRVs in 
CFO may enhance its exploration while preserving 
determinism so that the algorithm always yields the same 
results on runs with the same setup while effectively 
exploring DS. 


C. CFO Applications and Extensions 


CFO has been used to solve a wide range of practical 
problems, representative examples being training neural 
networks [2]; optimizing microstrip patch antennas [3,4]; 
designing broadband absorbers [5]; designing drinking water 
networks [6]; optimizing E-shaped laptop patch antennas 
[7,8]; designing tri-band slotted bowties [9]; analyzing 
financial systems [10]; assessing power grid reliability [11]; 
developing iris recognition software [12]; designing 
multilayer electromagnetic absorbers [13]; designing finite 
impulse response filters [14]; optimizing adaptive beam- 
forming antennas [15, 16]; designing UAV 3D flight paths 
[17]; designing reconfigurable [18] and tri-band handheld 
[19] RFID antennas; optimizing bandpass filters [20]; 
optimizing electromagnetic absorbers [21]; synthesizing 
antenna arrays [22,23]; increasing antenna impedance 
bandwidth [24,25,26]; investigating different objective 
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functions for broadband HF monopoles [27]; solving 
nonlinear circuits [28]; optimizing "smart" antennas [29] and 
a high voltage air-break disconnector [30]; improving two- 
level clustering for multi-view data analysis [31]; and 
optimizing transmitter locations for indoor optical wireless 
networks [32]. 


Proofs of convergence for CFO and extensions have been 
developed [33,34]. The algorithm has been implemented on 
a GPU using various neighborhood topologies [35,36]. 
Several CFO extensions and hybridizations have been 
developed, for example, adaptive DS and variable initial 
probes [37]; neural network ensembles using copula function 
theory [38]; real-time adaptation [39]; variable population 
size [40]; multistart CFO [41];  clustering-simplex 
hybridization [42]; CFO with differential evolution operator 
mutation [43]; adaptation based on stability analysis [44]; 
and hybridization with back-propagation simplex search 
[45]. CFO is included in several surveys of established 
optimization metaheuristics [46-51]. 


HI. WHY DETERMINISM MATTERS 


A distinction has been drawn between benchmark and 
real-world problems because real-world problems usually 
require a two-step solution whereas benchmark problems 
require only one. Optimization algorithms almost always are 
evaluated against benchmark suites for the very reason that 
the solutions are known in advance. An algorithm's 
effectiveness is measured by the quality of its solutions (how 
close it comes to the known answers) and by the required 
computational effort (usually the number of objective 
function evaluations). At no point is the fitness function 
itself an unknown, but this almost always is the case in real- 
world problems. The user must define a suitable objective 
function before any optimization can be done, and this added 
requirement makes solving the problem much harder because 
not all objective functions are well-behaved or for that matter 
yield good results even if they are. This idea is illustrated 
with some examples. 


A. A Typical Benchmark: Schwefel's Problem 2.26 


Schwefel's Problem 2.26 is a well-known benchmark used 
in many test suites. In an N} -dimensional decision space, 








Na 
this objective function is defined as f(x) => [x sinG/|x,))]> 
2 lrisin(y 


— 500 < x, < 500. It is highly multimodal with a single global 


maximum of 418.9829xN, at [420.9687]"* ([63] @ 
p.4467). The 2D maximum, which can be visualized, is 
837.9658 at (420.9687,420.9687) , but it is the 30D maximum 
of 12,569.487 that most often is used as a benchmark. The 
2D Schwefel's landscape, Fig. 1, is extremely multimodal 
which shows why this function is an especially challenging 
benchmark. 


Topology aside, testing against the Schwefel is 
straightforward because the problem is fully specified, both 
the functional form and its known solution. The only 
question is how well a particular algorithm performs against 
it. A stochastic algorithm might be run tens, hundreds, 
possibly thousands of times to determine how close it comes 
to the 30D Schwefel's known maximum of 12,569.487. 
Does knowing that an algorithm works well against this 
benchmark help when it comes to solving a problem where 
the objective function itself is unknown? Maybe, but 
probably not. 





Fig. 1. 2D Landscape of Schwefel's Problem 2.26. 


When the very form of the objective function is unknown 
the entire optimization process must be repeated for every 
candidate function. This is not a fundamental limitation, of 
course, but in many real-world cases as a practical matter it 
turns out to be one, and it cannot be overcome because as a 
practical matter computation times simply become too long. 
Evaluating the Schwefel is straightforward and quick 
because it requires only built-in functions (sine, square root, 
absolute value), but real problems usually are quite different. 


B. A Real-World Problem: Vehicle-Mounted Antenna 


Real-world problems often require an external "modeling 
engine” that frequently is far more computation intensive 
than using built-in functions. In antenna D&O, for example, 
Finite Element and Method of Moments codes are popular 
with their attendant long run times. This limitation renders 
statistical evaluations completely impractical. As an 
example, the run time for the 7,303-segment vehicle- 
mounted antenna model shown in Fig. 2 using the Numerical 
Electromagnetics Code (NEC4/MP, [66]) is 369 seconds on 
a 2.30GHz 8-core Core 17-3610QM machine running a 
multithreaded version of NEC4 under 64-bit Windows 7. An 
optimization run requiring, say, 1,000 evaluations would take 
more than 102 hours, well over 4 days! And this must be 
repeated for each candidate fitness function! 


Knowing that an algorithm is effective against Schwefel 
2.26 tells the practitioner absolutely nothing about how to 
formulate useful fitness functions and whether or not the 
algorithm will work well against those functions, whatever 
form they may take. Stochastic algorithms make answering 
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this and related questions very difficult and consequently are 
not well-suited to solving the real-world problems routinely 
encountered by practicing engineers. 


IV. A REAL-WORLD PROBLEM: MONOPOLE ANTENNA 


Because most real-world problems do not come with 
predefined objective functions, the first mandatory step is 
formulating one. How can the engineer know that a 
particular mathematical form is suitable to achieve the 
desired goals? How can the practitioner differentiate one 
function from another? These questions cannot be answered 
using a stochastic optimizer because every run yields a 
different result. 





Fig. 2. NEC Model of Vehicle-mounted Antenna. 


There simply is no way of knowing why the next run's 
results are different from the previous one's. Is it because the 
objective function is different? Or is it because of the 
algorithm's inherent randomness? Trying to answer these 
questions requires statistics generated over many runs, and 
even then the answer may not be correct because it is 
probabilistic. Deterministic optimizers like CFO avoid this 
conundrum entirely. If the objective function is changed, for 
example by adding a term or modifying one, then that 
change alone accounts for different results in subsequent 
CFO runs made with the same setup. Only deterministic 
optimizers allow different objective functions to be 
compared quickly and without ambiguity. 


A simple real-world problem illustrates how difficult 
defining a suitable objective can be and why using a 
deterministic optimizer like CFO can make a big difference. 
The example is optimizing a simple resistively loaded base- 
fed high-frequency (HF) monopole antenna on a perfectly 


electrically conducting (PEC) ground plane. The antenna is 
shown in Fig. 3 (see [27] for details). This problem is 2D 
with decision variables R and H (loading resistor value and 
height above the ground plane, respectively). The 
optimization objectives are maximum gain and as flat as 
possible an impedance bandwidth from 5 to 30 MHz. The 
relevant performance parameters are radiation efficiency, €; 
maximum gain, Gx? input impedance, Lin = Rin + JX in 
and voltage standing wave ratio relative to a reference 
impedance Zy)> VSWR (Z,) (the industry standard is 











Z =50Q). Three different objective functions are 
considered: 
AR Ss ae) (1) 
ow AVSWR (Z,) 
min (£) (2) 
R,H,Z,)= 
SaR, H, Za) |Z, — max (R, )|:|max (X „) 
fAR.H,Z,) = min (£) (3) 
woe" |Z, — max (R, )|: AVSWR(Z,)- [max (Xin) — min (Xin)] 





where AVSWR(Z,) is the difference between maximum and 
minimum VSWR over the 5-30 MHz HF band. 


At first blush it appears that each of these functions should 
achieve the desired objectives, but closer examination of the 
visualized functions shows this is not at all the case. 
Landscapes are shown in Figs. 4-6, and a cursory inspection 
reveals potentially serious problems. Function f, has a 
diffuse maximum close to the DS boundary. Its width may 
make locating the actual peak difficult, and its location near 
the boundary may inhibit effective exploration. f, is so 
pathologically spiky that many algorithms probably cannot 
efficiently locate its maximum. Only f, is well-behaved 


with a reasonably sharp maximum not too close to the DS 
boundary. 





PEC Ground Plane 


Fig. 3. Base-fed Loaded Monopole Antenna. 
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How can the antenna engineer decide which of these three 
ostensibly reasonable objective functions is best? In the 2D 
case the answer is simple: plot the landscape, look for the 
maximum. But 2D problems are rarely if ever the case and 
in any event do not require multi-dimensional search and 
optimization algorithms. What happens if the loaded 
monopole includes reactive loading in addition to resistance? 
Adding an inductor or capacitor increases the problem's 
dimensionality from 2D to 4D, the new decision variables 
being the value and location of the reactive element. Unlike 
the 2D problem the 4D problem cannot be visualized. The 
only way to decide which of many candidate objective 
functions is best is to make optimization runs with each one, 
a process that mandates the use of a deterministic optimizer 
like CFO in order to avoid the plethora of runs required by a 
stochastic optimizer. This example highlights the efficacy of 
algorithms like CFO when it comes to dealing with real- 
world engineering problems as opposed to benchmarks. 


V. ANOTHER REAL-WORLD ANTENNA PROBLEM: 
WIRELESS POWER TRANSFER 


Another real-world example involves optimizing a multi- 
arm helical antenna for wireless power transfer. This is a 4D 
problem with decision variables (i) helix length, (ii) 
diameter, (iii) number of turns, and (iv) number of arms. For 
the practitioner the first question is, How to combine these 
parameters? (Choosing a suitable fitness function requires 
investigating a variety of forms, which as discussed above is 
best done with a deterministic optimizer. 


The WPT optimization objectives are minimum VSWR, 
maximum power gain, and minimum helix volume at 200 
MHz. CFO was used as the optimizer, and after much cut- 
and-try the following less than obvious fitness function was 
determined to be a good one for achieving the objectives: 

1 
F= X (4) 
ka -|AGT -1|+|Z, - R, + A | 








Height (meters) 


800 2 
Fig. 4. Landscape of Objective Function f,(R, 0,50) - 


Note that VSWR does not appear explicitly whereas 
components of the antenna input impedance do, and that the 
value of Z, was determined by CFO using Variable Zp 


technology [64] (VZo) instead of being specified in advance. 
The other parameters are the free space wavenumber, k ; the 
largest dimension fully enclosing the helix, a; and NEC's 
free-space Average Gain Test, AGT (a measure of the 
antenna model's fidelity). The Numerical Electromagnetics 
Code (NEC) [65], is a widely used external Method-of- 
Moments modeling engine considered to be the "gold 
standard" for wire antenna modeling. 








400 


Load Resistor (ohms) = a “4 


800 `2 


Height (meters 


Fig. 6. Landscape of Objective Function f,(R, H,50) - 


Evolution of CFO's fitness is plotted in Fig. 7. The 
monotonic increase is typical of CFO's behavior even 
without elitism (preserving the best solution step to step) 
which was not included here. D,,., plotted in Fig. 8, is the 
average distance in DS between the best probe and all other 
probes. A value of zero means that all probes have coalesced 
on a single point. It is quite interesting that oscillatory 
behavior in Dayz apparently correlates with local trapping and 
perhaps can be used to signal that condition [54,55] and 
improve DS exploration. 


The optimized 200 MHz WPT helix is shown in Fig. 9. It 
has four upper and lower 1.537-turn arms wound with #12 
AWG wire. Its dimensions are 0.148m long by 0.089m 
diameter with the arms connected by 0.01m jumpers. The 
VZo optimized impedance is Zp = 95.670. VSWR//Zp is 
plotted in Fig. 10 (standing wave ratio relative to Zp), and the 
corresponding input impedance appears in Fig. 11. The 3D 
radiation pattern is shown in Fig. 12. Note that these data 
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assume PEC wires because the radiation efficiency with 
annealed Cu conductors is € = 99%. 


The optimized helix is electrically small (maximum 
dimension less than 0.1 wavelength at 200 MHz). Its VSWR 
is nearly perfect at 1.013:1. And the power gain is excellent 
at 1.88 dBi. The optimization objectives have been fully met 
because a suitable fitness function was employed whose 
form was determined by comparing candidate functions 
using deterministic CFO. Accomplishing this with a 
stochastic optimizer might require hundreds of runs for each 
candidate fitness function, which is clearly prohibitive for 
most real antenna D&O problems. 





Fitness (200 Mhz Multi-Arm Folded Helix) 
T T T T 


14 


Fitness 











Step # 


Fig. 7. Fitness Evolution, 200 MHz Helix. 


VI. BASIC CFO THEORY 


Basic CFO theory is presented in this section. Proofs of 
convergence and extensions developed by other researchers 
and are not discussed. 


A. Problem Statement 


Central Force Optimization searches a bounded N,- 
dimensional hyperspace Q for the global maxima of a 
fX, X253 Xy) 
determine max{ f(X): Xe QcR”}, f:QcR”—R with 
Q= {X x" <x <x, i=l, Ng} bounded 


region of feasible (allowable) solutions (decision space); 
f(X) a pointer to the objective function; and L 


=QuU f(X),X €Q. the problem's landscape. The x, are 


function Stated more precisely, 


being the 


the decision variables for each coordinate number i. Fitness 
refers to the value of the objective function f(x) at sample 


point X= (Xi X2- Xy,) in Q. There is no a priori 
information about the objective function’s extrema, that is, 


f (x)'s landscape (topology) is completely unknown. 


CFO searches © by flying "probes" through the space at 
discrete "time" steps (iterations). ach probe’s location is 





specified by its position vector computed from two equations 
of motion that analogize their real-world counterparts for 
material objects moving through physical space under the 
influence of gravity without energy dissipation. 
=> Na ‘ 
The position vector at step j is R? = De ees where 
k=l 


the x? are the probe's coordinates and ê, is the unit vector 

pae : er 5 
along the x,-axis. Indices p, 1< PSN, and j, 
O<JjJ<N, are the probe number and iteration number, 


respectively, where N p and N, are the corresponding total 


number of probes and total number of iterations (time steps). 





Davg 














Step# 





Fig. 9. 200 MHz 4-Arm Helix. 
(axis 0.05m, L=0.148m, D=0.089 m, #12AWG) 
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SWR (95.67 ohm) 200 MHz 
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Fig. 10. 200 MHz Helix VSWR // 95.67Q [PEC wires]. 
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Fig. 11. 200 MHz Helix Input Impedance [PEC wires]. 





y= 
Fig. 12. 200 MHz Helix Radiation Pattern [PEC wires]. 


B. Equations of Motion 


In metaphorical CFO-space each of the N p Probes is 


accelerated by the gravitational pull of the surrounding 
"masses." At step j—1 probe p experiences an acceleration 


given by the acceleration equation 


(Ri -RD 65) 


Np 
amO XUM ME) Mia) Ep 
j-l ja 


k#p 


which is the first of CFO’s two equations of motion. 
M? = fP, api. axei) is the fitness at probe p’s 
location at time step j—1. Every other probe at that step 
with it the fitness 
G is CFO’s gravitational 


Unit 


(iteration) has associated 
Mj ,,k=1,...,p—1, p+l,...,N,- 


constant, and U(-) is the 
1, z20 

U(z)= : 

(2) to See 


The acceleration a’, causes probe p to move from 


Step function, 


position RP, at step j—1 to position R? at step j 


according to the trajectory equation 


= > l . 

R? =R? +ã? At, j21 (6) 
J J 2 J 

which is CFO’s second equation of motion. Af is the "time 

interval" between steps during which the acceleration is 

constant. 


In the original CFO paper the counterpart of equation (6) 
above, that is, eq. (5) in [1], includes a velocity term that 
intentionally is excluded here. That term was set to zero in 
[1] as a matter of convenience because in the case of 
rectilinear motion it simply is an additive constant. But what 
eventually became clear is that including the velocity term is 
a mistake. It never should have been included in equation 
(5) in [1] in the first place because generally a probe's motion 
is not rectilinear. 


Trajectories instead are curvilinear which results in the 
acceleration and velocity vectors being in different 
directions. As an example, in the case of purely circular 
motion in real space the velocity vector is tangent to the 
circle while the acceleration vector is along the radius to the 
circle's center. The acceleration and velocity vectors actually 
are mutually perpendicular. This limiting case illustrates 
why, as a general proposition, the velocity term in Newton's 
real-space kinematic equations should not be included in the 
metaphorical CFO-space analogues - doing so erroneously 
changes the direction of a probe's acceleration. 


While CFO’s terminology mimics that of gravitational 
kinematics, it is important to remember that CFO, like all 
metaheuristics, is a metaphor, not a precise model of the 
physical system inspiring it. CFO's terminology has no 
significance beyond reflecting its kinematic roots, as, for 
example, does the factor /2 in equation (6). The gravitational 
constant, G, and time increment, At, have direct analogues 
in Newton’s equations in physical space, but the exponents 
@ and J do not. In keeping with CFO's metaphorical 


nature, these exponents are included to provide the algorithm 
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designer with added flexibility, say, by changing how gravity 
varies with distance or with mass or both if doing so creates 
a more effective algorithm metaphorical CFO-space. 


C. Mass in CFO-Space 


The concept of mass in CFO-space is very important, and 
it is quite different than it is in real space. In the Universe 
mass is an inherent, immutable property of matter. In CFO- 
space in stark contrast "mass" is a user-defined positive- 
definite function of the objective function’s fitness (not 
necessarily the fitness itself). 


In the acceleration equation (5) above mass by definition is 
MASS cq =U(M'_,—-M?,)-(Mi_,-M”.,)” , that is, the 
difference in fitness values raised to the @ power multiplied 
by the Unit Step. Of course, the algorithm designer is free to 
choose any functional form that works well. In this case the 
Unit Step is included to preclude negative mass because 
without it CFO's mass could be negative depending on which 
fitness is greater. 


Real mass in the physical Universe is always positive and 
consequently the force of gravity always attractive. But, 
depending on how it is defined, mass in metaphorical CFO- 
space can be positive or negative, which can lead to some 
very undesirable effects. Negative mass creates a repulsive 
gravitational force that flies probes away from maxima 
instead of toward them, thereby defeating the very purpose 
of the algorithm. While any definition of mass is permissible 
in CFO-space, only positive-definite ones allow the 
algorithm to perform as it should. 


D. Errant Probes in CFO-Space 


Another important consideration is the possibility that a 
probe’s computed acceleration may be so large that it flies 
outside Q. If any coordinate in the trajectory update 
equation (6) satisfies x; < x i or x% >x 


mx then probe is 
l 
errant because it enters a region of unfeasible solutions (ones 


that are not valid for the problem at hand). 


Errant probes arise in many algorithms, and the question is 
what to do with them. While there is a multitude of 
approaches, a simple empirically determined one is 
suggested here. On a coordinate-by-coordinate basis, probes 
flying outside DS are placed a fraction 0 < AF, <F „<1 


rep 
of the distance between the probe’s starting coordinate and 
the corresponding boundary coordinate. 

Fep is the variable repositioning factor discussed in more 
detail in [57-59]. Its value can be assigned by calculation, 
for example by incrementing it by a fixed amount A F , Or 
by specification, for example, using a fractions [62]. Not 
only does this scheme guaranty that all probes remain within 
Q , it probably improves CFO's exploration because F's 


rep 


value is uncorrelated with the landscape thereby adding a 
measure of pseudo randomness. 


VII. CFO PSEUDOCODE & IPD SPECIFICATION 


Central Force Optimization is simple and easily 
implemented. It comprises only the two equations of 
motion, (5) and (6) above, and may be implemented with the 
following steps: 


(1) (i) Compute (or assign) the initial probe distribution 
(IPD); (ii) Compute corresponding objective function 
fitnesses; and (iii) Assign initial probe accelerations (usually 
zero). 


(2) Probe-by-probe successively compute a new position 
in Q at the current time step using the previously computed 
(or assigned) accelerations. 


(3) Verify that each probe remains inside Q and make 
corrections as required. 


(4) Update the objective function fitnesses at each new 
probe position. 


(5) Compute accelerations for the next time step using 
the new probe distribution. 


(6) Loop over time steps. 


There are many ways to generate an IPD. The details are 
entirely up to the algorithm designer, and some IPDs will be 
better than others. It also is likely that different IPDs will 
work better with different classes of objective function 
landscapes. The "probe line" IPD [37] that seems to work 
well across a wide range of objective functions is described 
here. 


Fig. 12 is a schematic representation of a typical 2D 
IPD using probe lines. It comprises an orthogonal array of 


N 7 probes per dimension uniformly distributed on lines 
Ny 


parallel to the coordinate axes that intersect at a point along 
Q’s principal diagonal. The figure shows nine probes on 


each of two lines parallel to the x, and x, axes. The probe 


lines intersect at a point D=X,,.+7(X,,. —Xmi,) Where 


Na Na 
Y = mins and y = maxs are the diagonal’s 
Xaa e; X a= e; 8 

i=l i=l 


endpoint vectors. Parameter 0<y<1 determines where 


along Q's diagonal the probe line array is placed. 


Typical uniform 2D and 3D IPD's using variable probe 
lines (different Y values) appear in Figs. 13 and 14. It is not 


necessary that each probe line contain the same number of 
probes, and under some conditions different numbers may be 
more useful. As examples, equal probe spacing might be 
desired in a decision space with unequal boundaries or 
possibly to avoid overlapping probes. While the probe line 
IPD is completely deterministic, a random IPD could be used 
instead. a fractions, for example, could be used to create a 
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pseudo random IPD or they could be used to pseudo 
randomly place probes on probe lines. 


x 





yo" 
gor 


A. 
Np’Nd Probes per Axis 


we 


——f- 


Decision Space 
Fig. 12. Typical 2D Probe Line IPD. 


Representative CFO pseudocode using the variable probe 
line IPD appears in Fig. 15. Of course many adaptations are 
possible, for example, adding other loops to vary, say, the 
number of probes, the gravitational constant, or the CFO 
exponents. Other adaptations could include changing 
parameters or shrinking DS during a run in response to 
metrics such as rate of change or spread of fitnesses or 
perhaps probe convergence as measured by Davs. The 
possibilities are endless, and the efficacy of a particular CFO 
implementation is easily investigated because CFO is 
deterministic even with pseudo randomness. 


VIII. CONCLUSION 


Central Force Optimization has become an established 
metaheuristic for global search and optimization. It is 
particularly useful for real problems that require the 
definition of a suitable objective function before any D&O 
can be done. Stochastic algorithms make it difficult if not 
impossible to assess the utility of different fitness functions 
in real problems because each run yields different results, 
and there is no way of attributing the differences to the 
objective function being used or to the algorithm's inherent 
randomness. CFO eliminates this problem because it is 
completely deterministic, always yielding precisely the same 
result for the same run setup. CFO retains this characteristic 
even when it is implemented pseudo randomly using 
parameters whose values are precisely known but 
uncorrelated with the optimization problem's landscape. 
Adding pseudo randomness improves CFO's ability to 
explore a decision space, and some methods to that end have 
been discussed. CFO has been applied to a wide range of 
real world problems and has been improved and extended by 
many researchers. A bibliography of CFO-related work is 
included that may assist practitioners and researchers who 
are new to the algorithm. 





LOCATIONS OF 24 INITIAL PROBES FOR GP FUNCTION 
[gamma = 0] 











x1 
LOCATIONS OF 24 INITIAL PROBES FOR GP FUNCTION 
[gamma = .4] 











LOCATIONS OF 24 INITIAL PROBES FOR GP FUNCTION 
[gamma = .9] 








x1 


Fig. 13. Sample 2D 12-probes/probe line IPD’s. 
(y values: 0.0, 0.4, 0.9 top to bottom) 
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Fig. 14. Typical 3D 6-probes/probe line IPD’s. 
( y values: 0.0, 0.5, 0.8 top to bottom) 


For Y= Yoan t0 Yoop DY AY 

(a.1) Compute initial probe distribution. 
(a.2) Compute initial fitness matrix. 
(a.3) Assign initial probe accelerations. 
(a.4) Set initial Fai 


For j =0 to N, (or earlier termination criterion) 


(b) Compute probe position vectors 


R?,1< pS, leq]. 


(c) Retrieve errant probes (1 < DSN > ): 

If RPE, SP RPE =X + Fe (R a x") 

If RP 6, 2 xP RPE = xP + F(a R a ê) 

(d) Compute fitness matrix for current probe distribution. 

(e) Compute accelerations using current probe distribution 
and fitnesses [eq. (5)]. 


(f) Assign Ey g 





Fig. 15. CFO Pseudocode with Variable Probe Line IPDs. 
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